X-ray edge singularity of bilayer graphene 
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The X-ray edge singularity of bilayer graphene is studied by generalizing the path integral ap- 
proach based on local action which was employed for monolayer graphene. In sharp contrast to the 
case of monolayer graphene, the bilayer graphene is found to exhibit the edge singularity even at 
half-filling and its characteristics are determined by interlayer coupling. At finite bias the singular 
behaviors sensitively depend on the relative magnitude of fermi energy and applied bias, which is 
due to the peculiar shape of energy band at finite bias. 
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Introduction- Physics of graphene systems is very 
richi A monolayer graphene has massless Dirac fermion 
band whose density of states at Fermi energy vanishes 
at half-filling, and this is responsible for many of its 
semimctal properties X A bilayer graphene is composed 
of two monolayer graphenes connected by interlayer cou- 
pling in Bernal stacking structured— The band struc- 
ture of a bilayer graphene can be controlled externally 
by applying a gate bias between two layers^ The (low 
energy) band of unbiased bilayer graphene is that of mas- 
sive Dirac fermion 2 without gap between electron and 
hole bands, while a bias opens energy gap between them. 

Meanwhile the transport properties of graphene sys- 
tems are intensively studied, the optical studies also re- 
vealed many interesting physical properties. e Among di- 
verse optical probes we focus on the (near-edge) X-ray 
absorption spectroscopy^ and X-ray photoemission. The 
physics of X-ray edge problem of fermi liquids is very 
fascinating and well understood4~~ The incident high- 
energy X-ray photon excites a deep core electron, leav- 
ing behind a (positively charged) core hole which can be 
treated to be immobile in many cases.^ If the excited 
electron escapes into vacuum this is called X-ray photoe- 
mission process, and if the electron cannot escape this is 
referred to as X-ray absorption process. Now the fermi 
sea of conduction electrons forming fermi liquid reacts to 
this suddenly created (namely, time-dependent) potential 
which can be assumed to be well localized for the energy 
in the vicinity of absorption edge^ The conduction elec- 
trons interact with the potential in two distinctive ways: 
the excitonic processes 1 ^ which is essentially attraction 
between hole and conduction electrons and the orthog- 
onality catastrophe^ which means a vanishing overlap 
between the ground state wavefunctions of conduction 
electrons before and after the creation of the hole. Both 
excitonic processes and orthogonality catastrophe are sin- 
gular near X-ray absorption edge for fermi liquids in the 
sense that there are divergences in the perturbative ex- 
pansion, so they require nonperturbative treatments-that 
is solving a particular type of singular integral equation 
exactly.-^ This singular feature is strongly dependent on 
the fermi liquid nature of conduction electrons, especially 
the finite density of states at fermi energy. Then it is nat- 



ural to ask what will happen for a graphene which has a 
semimetal character at half-filling. 

Motivated by the above observation the author has 
studied the X-ray edge problem of monolayer graphene 
in ref.[3] (hereafter called (I), and see also ref.[l5[)- The 
results for the monolayer graphene can be summarized 
as:(l) at half- filling the edge singularity is eliminated due 
to the vanishing density of states at fermi energy (2) away 
from half-filling the edge singularity revives, and the ex- 
ponent of power law behavior of the singularity depends 
on fermi energy non-monotonously. 

In this Brief Report we report the results of on the X- 
ray edge problem of bilayer graphene which are obtained 
by generalizing the approach of (I) , highlighting the dif- 
ferences with the case of monolayer graphene. The most 
important finding is that the X-ray edge singularity per- 
sists for the bilayer graphene even at half-filling, and this 
is essentially due to the nontrivial matrix element of 1- 
particle Green's function which stems from the interlayer 
coupling. Also for the case with finite bias the edge sin- 
gularity shows up, but its characteristics very sensitively 
depends on the relative magnitude of fermi energy and 
applied bias. This feature is due to the peculiar shape of 
the energy band with finite bias. Since the mathematical 
framework of approach in this report is almost identical 
with that of (I) , we avoid duplicating the essentially iden- 
tical mathematical expressions, instead we refer readers 
to (I) for more details. 

Setup- The honeycomb lattice of monolayer graphene 
is spanned by the following basis vectors: ai = £(1, 0) 
and a 2 = where £ = 2.46A is the lat- 

tice spacing. Each site of A sublattice is connected 
with three B sublattice sites (assuming a A site is at 
(0,0)): n = £(0,^-), r 2 = i{-\,-JL=), and r 3 = 

£{\,— )• A bilayer graphene has Bernal stacking 

structure^, such that A site of top layer(=layer t) is sit- 
ting directly on the top of B site of bottom layer (=layer 
b). The unit cell of bilayer graphene contains four sites, 
and we choose the ordering of the sites according to 
(bA, bB, tA, tB) = (1, 2, 3, 4)M The tight-binding Hamil- 
tonian in Slonczewski-Weiss-McClure parameterization 4 
includes one intra-layer hopping 70 and three inter-layer 
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hoppings 71,73,74. The most important characteristics 
of energy bands are determined by 70 and 71^ hence 73 
and 74 will be ignored in this paper. The values of 70 
and 71 are 70 ~ 3.2 eV, 71 ~ 0.4 eV£ We also take into 
account the difference in site energy u between two lay- 
ers which can be controlled externally by applying bias 
voltage, 5 Then the tight-binding Hamiltonian of bilayer 
graphene in momentum space takes the following form 
(spin indices are suppressed): 



Ha 



£4 

k 



■u/2 7o A k ■ 

70A* u/2 71 

71 -u/2 70 Ak 

. 7o A* -u/2. 



Ck, (1) 



where Ck = (ci,k, C2, k , C3k, C4,k) 4 (t is matrix transpose) 
and c^k is the electron destruction operator with layer 
and sublattice index i. Ak = 53i=i e lk ri . The diagonal- 
ization of eq.([T]) gives four energy bands ±_B±k, where 
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T + W 



(2) 

Note that E + ^ > 71. ±E + ^ is the high energy electron 
(hole) band (often called dimer bands). X-ray edge sin- 
gularity is low energy process involving conduction elec- 
trons near fermi energy. Therefore, we will be mostly 
interested in the low energy bands ±-E_ )k which is gap- 
less when u — 0. To see this, note that Ak vanishes 
at six corners of Brillouin zone, and among them only 
two are distinct: K + = (|f,0) and K_ = K + . This 
naturally introduces the valley index ±. Expanding 
Ak around each valley K±, we obtain (for small q x ,y) 



A 



K++q 



2 v 



Qx+iq y ) and A K _+ q 



iq v ) 



Now in the small q (continuum) limit, the Hamiltonian 
(fTJ) becomes 8x8 matrix ( 16x16 if spin included) 



H 



f d 2 q 


'* f + (q)" 




k + " 




*+(q)" 


1 (2tt) 2 


* f -(q). 




a:_ 




tt_(q)_ 



(3) 



where \&(q) = v VC^ (V is the total volume of lattice) 

and *±(q) = [^±i(q), V'±2(q), V>±3(q), ^±4(q)]*- K + is 
the following 4x4 matrix for the valley K + : 




K, = 



^"fal/h and n(q) 



(4) 



v(—q x +iq y ). The matrix 
can be obtained through AT_ = 



where v 

K_ for the valley K_ 
K+{q x -> -q x ). 

In this report we neglect high energy processes involv- 
ing the dimer bands ±i? + q . We will take the 71 ~ 0.4eV 
to be the highest energy scale of our problem, and the low 
momentum region is defined by the condition vq <C 71 

(q = \/Qx + 1y)- It is also fairly reasonable to assume 



that 71 3> \u\. In realistic experimental conditions, the 
fermi energy fi is also much smaller than 7! . If u = then 
EL.q ~ v 2 q 2 /ji in the low momentum region, which is 
the energy dispersion for the massive Dirac fermions. If 
u 7^ the low energy electron band has a local maxi- 
mum at q — with i?_, q= o = and a global minimum 

_ M 



at q nln — , u y -^f +u 



= &JW£ withi?_ 



71 



Since 



1 2 v^F - 

M ^ 71 , the above implies that the low energy band 
is almost flat in low momentum region. For conduction 
electrons to be available (electron doping assumed) with 
finite u, the fermi energy should be highter than EL. 9min . 
We will assume fi > \u\/2. 

The Hamiltonian of the (almost immobile) deep core 
electron is taken to be Edd^d, where Ed < is the core 
level energy (d is the destruction operator of deep core 
electron). The scattering potential by deep core hole is 
assumed to very locals, so that it is diagonal in layer 
and sublattice indices. Then the interaction Hamiltonian 
between conduction electrons and deep core electron is 
given by 



H u 



d 2 x¥(x)6(x)V^(x)dU, 



(5) 



where V is a 8x8 matrix representing intra- valley (Vo) 
and inter- valley (Vi) scattering^ 



V = V Ig + Via x <g> I 4 , 



(6) 



where a x is the 2x2 Pauli matrix acting on the valley 
space, and V\ is assumed to be real for simplicity. I4 and 
Is is the identity matrix acting on layer-sublattice space 
and valley-layer-sublattice space, respectivley. 

Local action approach- Since the scattering process oc- 
curs only at x = our problem is essentially local, and it 
is advantageous to formulate the problem in a local way 
by integrating out the bulk degrees of freedom except for 
the one at x = i 14 i 18 i 19 The resulting local action in 
imaginary time is (see (I) for more details) 

S[ V ,d}=-J2 I dTdr'n a {r)[G^]-J{T-T') Vh {T') 



drd{T)d{T)[fj(r)Vr]{T)] + / drd(d T - ur)d, (7) 



where a, b are the valley-layer-sublattice indices and lot = 
/i — Ed is the (bare) threshold energy. r\ a = ^ a (x = 0) 

is the local degrees of freedom. G a Q b \r — r') is the local 
Green's function (zero temperature assumed hereafter): 



-<*„(£ =0,r)*J(f =0,r')) 
ded 2 q„ e(7 



(2tt) 3 



G^feq), 



(8) 
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where 

.(0) 



(ie + n)U -K+ 

(ie + fi)h - K- 

G+(ie,q) 
G-(ie,q) 



(9) 



The matrices K± are defined in eq.Q and G ± is the in- 
verse matrix of (ie + fi)!^ — K±. Upon angle interal of 
q, it is straightforward to check J dqG + = J dqG~ . The 
most essential difference between monolayer and bilayer 
graphene lies in the eq.® as will be discussed in de- 
tail below. The explicit momentum integral shows that 
only G±(t) with (ij) = (11, 44, 22, 33, 23, 32) are non- 
vanishing and that Gj 3 = G^. 

Nozieres and De Dominicis (ND)'s solution— Let us 
consider a single species of fermion whose (unperturbed) 
local Green's function in the long time asymptotic limit 
is given by G(°'(t) = — p/r, where p is the density of 
states at fermi energy. Then the Green's function in the 
presence of time-dependent potential — V9(t\ — t)9(t — 
t 2 ) (0(x) is step function and V is a positive constant) 
satisfies the following singular integral equation, 

g& a = g (o) (c - a + f T1 dT G(0) « - t)(-v)g(t, o. 



The solution of eq. pO 
the long time limit is 



(10) 

which is asymptotically exact in 



Gnd(£,£'|ti,t 2 ) = [cos 2 6\GW(S- 



K-r 2 )(n-r) 



■r 2 ). 
(11) 



fl2) 



with the scattering phase shift given by 

5 = tan -1 [7rFp]. 

Correlation functions- In our case there arc many 
species of fermions, so that the eq. tfTU]) becomes a ma- 
trix equation: G (0) ->■ and V -> f [see eqs. (jrJI5|) ]. 

Also, the exact solution cannot be obtained for G^(r) 
of general form. 

The photoemission is related to core hole (Matsubara) 
Green's function [see eqs. (28,29) of (I)] 

D(t u t 2 ) = (T T d\ Tl )d(T 2 )) 



G (o)v) 



(13) 



= exp Tr ln(l 

where V is the time-dependent potential (matrix) 

V( T ) = -V6(n - t)9(t - r 2 ). (14) 

The X-ray absorption is related to the hole-conduction 
electron response function [see eqs. (38,39) of (I) ] 

F(n,Ti) = £ (T T dt(ri)»; (Ti)^(7a)d(7a)) 

a . b 

= -e(n-T2)D(n,T 2 )J2Gab(n-T 2 )(l5) 



where G a b(T~i — T2) is the solution of ea. (110[) in matrix 
form. Note that £ — > t\ — r c and £' — > t 2 + r c should 
be taken to obtain G a b in eq. (|15j) (t c is a short time cut- 
off). The core hole Green's function £>(ti,t 2 ) can be also 
obtained from G a b(T\ — r 2 ) by parameteric integral] 12 i 20 
Therefore, the task is (1) to compute the unperturbed lo- 
cal Green's functions of eq.©; (2) substitute them into 
the matrix version of eq. (I10[) and find its solution (if pos- 
sible); (3) do parameteric integral to find Z)(ti,t 2 ), and 
finally obtain F(ti,t 2 ). For the computation of unper- 
turbed local Green's functions of eq.© it is convenient 
to consider two cases, u = and u ^ 0, separately. 

Results for the unbiased case - For u = 0, G n = G44, 
G^ 9 = G^ 3 hold. In the long time limit we obtain 



Gf 1 (r) = 



GUt) 





(71 -m) 1 


1 




8ttv 2 t t 


87TD 2 T 2 




(71-M) 1 






87ra 2 T ~ 


8t™ 2 t 2 




r m 


1 






8ti"i; 2 t 2 




1 ^ 


l-2e<" 






87ti; 2 t 2 



T > 0, 
T < 0. 



T > 0, 
T < 0. 



(16) 



(17) 



G 23 (r) is identical with G 22 (t) up to sign. Eqs. (|16ll7[) 
should be compared with those of monolayer graphene. 



G„ 



,(t > 0) 



47Ti; 2 7 

M 

47T"U 2 7 



47rt; 2 r 2 



t > 0, 
t < 0. 



(18) 



It is clear G^r) is essentially identical with G m0 no(T) 
up to a factor of 2 (which comes from neglected dimer 
band of bilayer). However, G 1:l (t) is qualitatively dif- 
ferent from G mono (r). The difference is most marked 
at half-filling /i = 0. In this case, the leading term of 
G 11 (t) is proportional to 1/r, while that of G mono is 
proportional to 1/r 2 . It is easy to check that when the 
unperturbed local Green function is proportional to 1 /r 2 
there is no edge singularities in perturbative expansion 
[see eq. (46) of (I)]. Therefore, we conclude that the 
X-ray edge singularity persists for the bilayer graphene 
even at half- filling. Next let us compare Gj^t) with 
G^(r). The difference at half- filling jj, = is manifest, 
and we also note that even away from half-filling (/i 7^ 0) 
71 ^> fi holds, so that the most dominant contribution 
comes from the —ji/8ttv 2 t piece of Gf^r) (and G 44 (r)) 
which is the premise for the ND solution. This appear- 
ance of 71 is due to the non-trivial matrix element of 
1-electron Green's function (0\ipn(x = 0) | q) . Then the 
matrix integral equation is almost diagonal except for the 
potential matrix V. However, this can be dealt with the 
ansatz [ see eqs. (20,21) of (I), there t x should be corrected 
toI 2 ], 



G(r) 



h 



• G'(t)+<t x ®G"(t), 



(19) 



where G' and G" are the diagonal 4x4 matrices whose 
[22] and [33] entries are absent, and [11] and [44] entries 
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are identical. Then the solution of integral equation is 
given by [see eq.(25) of (I)] 



r'<" u^rW 1 ) cos 2 5+ 

^11,44^ > °) ~ 2 T 

1 (~7i)cos 2 5- 

2 T 



25+ /n 



2(5- /tt 



where the scattering phase shift S± is given by 



5± =tan- 1 [7r(y ±^i) 



7i 
8irv 2 



(20) 



(21) 



Note that is the counterpart of the density of states 
at fermi energy for fcrmi liquids. The core-hole Green's 
function is given by 



D(t > 0) - e 



1 



{t/t c ) n °( 5 M 2 



(22) 



where N c = 4 ( ±11, ±44 ) and S 2 = (5\ + S 2 _)/2. lu^ 
is the renormalizcd threshold energy. Finally the X-ray 
absorption intensity which can be obtained from F(t) by 
Fourier transform and analytic continuation is given by 
(E c is energy cutoff) 



I(lu) ~ 0(uj — Wy) 



7i 

inv 2 



OJ — UJj 
E r . 



-2<5 + /7r+AT c (<5/?0 



(23) 

Results for the case with finite bias- At finite u, asym- 
metry between top and bottom layer is expected, and 
this is reflected in the structure of local Green's func- 
tion. In this case the unperturbed local Green function 
is determined by the momentum region around fermi line. 
Recalling the assumption \x > u/2, we find that the fermi 
momentum is given by 



2 1 
% = — 



,i 2 + « 2 /4 + Ju 2 p? + 7l 2 (/i 2 - (u/2) 2 ) . (24) 



Then have E. 



v 2 (q 2 — go)-, where c is a dimension- 



less constant given by 

,2 




2 2 2 
U V Qq 



(25) 



which becomes very small in the limit \i — > u/2. This fea- 
ture originates from almost flat band strucutre at finite 



u which was mentioned previously. The long time limit 
of unperturbed Green's functions can be computed using 
the Laplace method of asymptotic analysis: 



GfAr) 



1 



c Wr u ' 



(26) 



where are dimensionless constants whose detailed 
forms do not concern us here. The constants just show 
that 33 23 can be ignored compared to the G^ 1Ai . 
When fj, is close to u/2, only Gfi gives a dominant con- 
tributions. In this case, the results eqs. (|22!23l) still apply 
with the modification 



7i 



fj,zu 



Nr. 



--11). 



(27) 



Owing to the factor c of eq. (|25[) , a very interesting varia- 
tion of correlation function as a function of fermi energy 
and bias is expected. 

Remarks and summary- The near edge X-ray absorp- 
tion experiment for graphene systems was reported in 
refs.[21]. The K edge singularity around 283 eV is clearly 
visible. However, the fermi energy and the bias depen- 
dence, which is our main point, was not studied in this 
experiment. For more precise comparison, we have to 
include the effects of the (ignored) dimer band, the con- 
tribution from C Is a* transition, the band dispersion of 
core level, and various broadening effects such as lifetime 
and temperature, and full treatment of these are beyond 
the scope of this Brief Report. 

In summary, we have studied the X-ray edge singu- 
larity of bilayer graphene by generalizing the approach 
employed for monolayer graphene. The bilayer graphene 
exhibits the edge singularity even at half-filling. Also 
at finite bias, the singular behaviors are found to depend 
very sensitively on the relative magnitude of fermi energy 
and applied bias. 
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